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ABSTRACT 

Context. The standard picture of accretion is a steady flow of matter from the disc onto the young star - a concept which 
assumes the star-disc system to be completely isolated. However, in a dense cluster environment star-disc systems do 
interact gravitationally. 

Aims. The aim here is to estimate the encounter-induced accretion rate in an ONC-like environment. 

Methods. Combining simulations of the cluster dynamics with simulations of the effect of encounters on star-disc systems 

we determine the likelihood and degree of encounter-triggered accretion processes. 

Results. We show that accretion bursts triggered by encounters of star-disc systems are common in young dense clusters 
like the ONC leading in the outburst phase to typical accretion rates of 10~ 7 -1CP 4 Mo/yr. Up to a third of stars presently 
in the Trapezium region accreted at least 1% of their disc mass via this mechanism in the last lMyr. Accretion of over 
6-7% of the disc material can occur in a single encounter. Despite losing their discs quickly, the total percentage of disc 
matter accreted per star is largest for the massive stars. 

Conclusions. Supplementing the steady accretion flow there exist episodic periods of high accretion in dense cluster 
environments. Due to their high accretion rate these processes should be observable even now in some of the low-mass 
stars in the ONC. 
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1. Introduction 

The theory of accretion has been extensively studied 
in the context of t h e standard a c cretion disc model 
(jPringle fc Rees 19721 : IShakura 1971 iPringle 19811 1. The 
main idea is that of gas moving inward due to radial 
turbulent transport of angular momentum to the outer 
disc regions. Currently the magneto- rotational instability 
(Bal bus fe Hawlev 19911 ) seems to be able to provide the 
required degree of turbulence in accretion discs. 

The disc evolution that follows is not completely un- 
derstood. In the later phases accretion happens exclu- 
sively from the disc material, so linking disc development 
and accretion directly. The often-used theoretical mod- 
els of viscous disc evo lution (|Lvnden-Bell fc Pringle 1971 : 
lHartmann et al. 1998f ) treat the disc as an isolated smooth 
axisymmetric structure that evolves due to an unspeci- 
fied source of temporally independent, turbulent viscos- 
ity. These standard models assume a possibly decreasing 
with time but relatively steady accretion process. Recently 
models were suggested where accr etion is caused by grav- 
itational instabili ties in the disc (jVorobvov fc Basu~ 2006: 
iBolev et al. 20061) leading to alternating phases of high and 
low accretion. However, for these instabilities to occur the 
disc has to be relatively massive (rrid > 0.1 M*), which 
means it will predominantly occur very early in star forma- 
tion. 

From observations the general picture emerges that 
the disc fraction is a strong function of age, decreas- 
ing from ~ 80% for clusters at 1 Myr down to few 
discs for cluster s older than 10 Myr (H aisch et al 2001: 
iHernandez 20071 ). The typical acc retion rate for a young 
star is ~ 10" 8 - lO- lo M /yr ([Hillenbrand et al. 1993 



iGullbring 19981 : lHaisch et al 20011 ). However, the devel- 
opment of the disc and the accretion rate of in- 
dividual sour ces in the same environm ent can differ 
considerably ( Fl aherty fc Muzerolle 20081 ). The accretion 
rate depends on the stellar age, mass and environ- 
ment. As expected from viscous disc models a de- 
crease in accretion rat e s with stellar age has been found 
(|Hartmann et al. 19981 iMuzerolle et al. 2000D . A matter 
of debate i s the dependence o n the stellar mass of the 
form M a (ICalvet et al. 2004 ISicilia-Aquilar et al. 20061 
iNatta et al. 2006D . possibly with a ~ 2. The effect of the en- 
vironment possibly manifests itself in th at Herbig Ae stars 
have on average higher accretions rates (|Lopes et al. 2 006) 
than low-mass stars in Taurus and Ophiuchus. 

Here we investigate an accretion process that will be 
present in dense cluster environments in addition to the 
steady accretion - the passage of a star inducing a simul- 
taneous transport of disc matter outwards and inwards via 
spiral arms, leading to accretion. We will show that the con- 
sequences are short bursts of high accretion (10~ 7 - 10~ 4 
M /yr), inevitably occuring in dense clusters and observ- 
able due to these high accretion rates. 

Although this process has properties in common 
with accretion induced by gravitation al instabilities 
(IVorobyov fc Basu 20061 iBolevet al. 2006f) and in wide bi- 
naries (jBonnell fc Bastien 1992t ). the cluster dependence is 
unique to the process treated here. 

For encounters to play an important role the cluster 
has to be dense. A typical example for such an environ- 
ment is the Ori on Nebula Cluster(ONC). In an earlier paper 
()Pfalzner 2 006). we demonstrated that encounters cause a 
3-5% specific angular momentum loss in the ONC rising to 
10-15% in the dense inner Trapezium region. Since specific 
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angular momentum loss is a prerequisite for accretion we 
suggested that in the final star formation stages an addi- 
tional growth mechanism for massive stars exist (cluster- 
assisted accretion) as massive stars lose their specific angu- 
lar momentum to a higher degree than low-mass stars. Here 
we study encounter-induced accretion directly by measur- 
ing the amount of matter kicked into the inner disc areas. 
Obviously this is not strictly accretion, but represents a 
often-used method in circumstances where computational 
expense does not allow to further resolv e the inner disc 
(jBate et al. 2002t IVorobvov k. Basu 20051 ). In the following 
the term "accreted" will be used in this sense but the lim- 
itations of this approach should be kept in mind. 



2. Method 

We combine simulations of the cluster dyna mics performed 
with the Nbody6+- 1- code (|Spurzem 19991 ) with tree code 
simulations to study the effect of encounters on star-disc 
systems as described in Pfalzner (2006). Here we determine 
the likelihood and degree of encounter-triggered accretion 
processes. For simplicity all stars of the cluster are assumed 
to be initially single, the effect of gas and the potential 
of the background molecular cloud OMC 1 arc neglected. 
Cluster models were set up with a spherical density distri- 
bution p(r) oc r~ 2 and a Maxwell-Boltzmann velocity dis- 
tribution. 9 1 C Ori was placed at the cluster centre and as- 
signed a mass of 50 M Q . For all other stars the masses were 
gen erated random ly according to the mass function given 
by (|Kroupa 2002D in a range 5OM > M* > O.O8M . The 
Nbody6+- 1- code determines th e trajectories of the stars to 
a high precision (|Aarseth 20031 see). Information of all per- 
turbing events of each stellar disc was recorded during the 
simulation, i.e. both masses, the relative velocity and the 
eccentricity. More details about the encounter simulations 
can be found in the appendix. The observed peak density 
of the Trapezium region of 4.7 xlO 4 stars/pc -3 was repro- 
duced in our simulations within 20 %, as was the stellar 
density distribution of the ONC. We chose the ONC to be 
in virial equilibrium. 

The code applied here cannot describe the disc dynam- 
ics simultaneously with the accretion process. The mass 
density - as the required number of simulation particles - 
would be too high close to the star and so would be the 
required temporal resolution. In any case the real accre- 
tion process will be rather complex. The main question is 
whether perturbations to the inner disc area can be induced 
by an encounter. Whether this is done directly by a parti- 
cle stream or only some related mechanism is of secondary 
importance. Therefore we resort to the approach of Bate 
et al (2002) and Vorobyov & Basu (2005), by computing 
the amount of matter reaching a sphere around the central 
star (here 1AU). We model only gravitational interactions, 
hydrodynamical effects have been neglected. In principle 
the latter could hinder matter reaching the inner disc ar- 
eas due to pressure gradient. However, here the velocity of 
the instreaming matter is so high that this is in most cases 
not a problem. Nevertheless the accretion values attributed 
here should only be regarded as relative values, determining- 
actual accretion rates would require a more sophisticated 
approach. The present study only aims at determining in 
which context encounters are likely to trigger accretion pro- 
cesses. 



3. Encounter results 

First we determine how much disc mass is accreted depend- 
ing on the encounter parameters. In principle, the star can 
not only accrete matter from its own disc but from the disc 
of the passing star as well. We find this only happens in 
very close or penetrating encounters and so this process is 
not included in the following. We performed a parameter 
study of coplanar, parabolic encounters, which shows that 
in a single fly- by up to 6-7% of the disc mass (Fig. IA.1I and 
Tabic lA.ip is accreted. In general much more mass is lost 
than accreted, with typical ratios of m. accretec ;/TO/ ost ~ 2-5 
%, which can rise to ~ 10% in some cases. Even fly-bys as 
distant as 10 ra, where r<j is the disc size (150 AU), can 
induce accretion provided that the mass ratio M^/Mf of 
the two interacting stars is high enough, where Mj* indi- 
cates the mass of the considered star and M£ the mass of 
the passing star. For Mj/M-j" < 1 the amount of accreted 
matter simply increases the closer the encounter, but for 
very small mass ratios, M^/M^ < 0.1, the accreted matter 
never exceeds 1% of the disc mass regardless of how close 
the star passes. For M^/Mj* > 1 the higher mass passing 
star increasingly captures matter that otherwise would have 
been accreted by the primary. Note, although the above re- 
sults were obtained for discs with masses ma = 0.001M Q , 
they are applicable for ~ md > 0.5Mj\ For a discussion see 
Appendix. 

Typical disc masses are rrid ~ 0.01M© and accretion 
bursts last normally ~ 10 2 - 10 3 yrs. So 5% accreted disc 
mass is roughly equivalent to an averaged accretion rate 
of 5 xl0~ 7 — 10~ 6 M Q /yr. However, there will be a wide 
range of induced accretion rates, as the disc mass depends 
strongly on the stellar age and mass. 

So far we considered only prograde encounters, but what 
happens in retrograde encounters? As long as the passing 
star does not penetrate the disc, no matter becomes ac- 
creted. However, encounters closer than 0.5rd can lead to 
higher accretion rates in retrograde than in prograde en- 
counters. As distant encounters are more common using the 
accretion rates from prograde encounters leads probably to 
a somewhat overestimated accretion. 

Cluster dynamics simulations show that especially in 
the dense inner regions star s often undergo a series of en- 
counters (jOlczak et al. 20061) . We simulated consecutive en- 
counters with the same parameters but using the resulting 
density distribution as input for the next encounter. Each 
encounter reduces the disc mass and alters the density dis- 
tribution. The result is that in subsequent encounters less 
matter is accreted. To approximate this, we use the results 
from single encounters but relative to the reduced disc mass 
determined using Eq. (4) from Olczak et al. (2006). 

In the following we combine abo ve encounter resul ts 
with cluster simulations described in (|Qlczak et al. 20 06V 

4. Cluster-induced accretion 

For ONC conditions nearly all cluster-induced accretion 
processes happen within 0.2pc around the cluster center. 
Here the high stellar density leads to a large number of 
encounters(> 500 within 0.03pc compared to < 50 outside 
0.2 pc). 

After ~ 0.1 Myr about a third of all stars presently 
in the dense Trapezium region have accreted at least 1% 
of their disc mass (see Fig. [I}. In other words, these stars 
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Fig. 1. Percentage of stars currently in Trapezium and 
ONC region with < 1% of the disc mass accreted as function 
of time and the same for stars that transversed these areas 
in their past. 
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Fig. 2. a) Mass accreted per encounter in the Trapezium 
region and the entire ONC averaged over the first 3 Myrs 
and b) mass accreted per star during the accretion phases 
shown as a function of primary mass in units of disc masses. 



have shown signs of very high accretion (> 10~ 7 M Q /yr) 
due to interacting with neighbouring stars in the past 1-2 
Myr. The stellar content fluctuates, so looking at stars that 
have visited the Trapezium region in the past, this fraction 
will increases to ~ 70% by 3 Myr. 

Fig. [2^i) shows the accreted matter per star as a func- 
tion of the primary mass Mj" during the first 3Myr. In the 
Trapezium region the accreted matter rises by a factor of 
~ 1.5 from low-mass to massive stars. In the entire ONC 
the dependency is much stronger as many more low-mass 



stars never undergo an encounter leading to accretion. Here 
the average accreted mass per star is 10-15 times larger for 
massive than for low-mass stars. 

We find that the accreted mass per accreting encounter 
is nearly constant up to M± ~ 5 M Q and rises for higher 
masses. The average number of accreting encounters is also 
higher for massive stars. In the Trapezium region a 20 M Q 
star has on average more than twice the number of accreting 
encounters than a 0.1 M Q star; this factor increases to 10 
times when considering the entire ONC. However, the most 
massive stars lose their discs quickly (typically ~ 5 xlO 5 
yrs) so that although they undergo encounters there is no 
matter left to be accreted. 

Nevertheless, the accreted mass is highest for the most 
massive stars. The mass dependence can be most clearly 
seen looking at the average accreted mass in the accretion 
phase - meaning the time until the star has completely lost 
its disc. Fig. [2Jd) shows that the accreted mass is more than 
30 times more for massive stars than low-mass stars. 

The number of accreting encounters decreases strongly 
with cluster age, so are there currently any observable 
encounter-triggered accretion bursts in the ONC happen- 
ing? Although each individual low-mass star is less likely 
to undergo an accreting encounter than massive stars due 
to their relative abundance, one is more likely to observe 
a low-mass star having such an encounter than stars with 
higher mass (see Fig. [3]). For the low-mass stars the number 
of encounters decreases from several hundred in the first 10 4 
yrs to 5-10 such events in an timeinterval of 10 4 yrs at the 
current age of ~ 1 Myr. As most encounter triggered accre- 
tion processes last 10 3 - 10 4 yrs, this means that one can 
expect to find currently between 1 and 10 low-mass stars in 
the ONC that show high accretion rates due to encounter 
processes. 

5. Discussion and Conclusion 

There is one observed systems that displays the link 
between encoun ter and accretion rate. The morphology 
and kinematics (|Cabrit et al. 200lf ) of RW Aurigae A and 
B strongly suggest tidal stripping of the primary disc 
by a recent fly-by that is occurring. This system dis- 
plays a high accretion rate of ~ 2-10 x !0~ 7 M©/yr 
(jBasri fc Bertout 1989t lHartigan et al. 1995T) while having 
a particularly low disc mass of m c d ~3 x 10 -4 M©. Using 
these parameters our simulations show the likely total ac- 
creted mass to be 0.3 - 4 % of the pre-encounter mass, 
equivalent to a maximum induced accretion rate of 1-10 
- 1.2 -10~ 5 M©/yr assuming a pre-encounter disc mass 



< 10 x m^) and that most mass is accreted within 



~ 10 years (although high accretion for several 100 years 
exist). This accretion rate is consistent with the observed 
accretion rate of ~ 2 — 10 • 10~ 7 M Q /yr. So encounter- 
induced accretion is probably occuring in this system. 

However, RW Aurigae is not located in a dense cluster, 
so either the rare case of an encounter in a low-density 
region or a strongly elliptical wide binary is being ob- 
served. Even in the latter case, the results described in 
Section 3 are still applicable as the effect of a parabolic 
encounter on the induced accretion is very similar to that 
in a strongly elliptical wide binary. Episodes of rapid ac- 
cretion due to tidal disturbances to the disc induced by a 
binary companion were suggested by Bonnell and Bastien 
(1992) in the context of FU Orionis flareups, which are 
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jet sources are frequentl y (in 85% of c ases) found to have 
close stellar companions (|Larson 20021) . 



6. Summary 

We have demonstrated that encounter-induced bursts of 
accretion are common phenomena in the centre of young 
dense clusters. Successive encounters can lead to over 10% 
of the disc mass being accreted typically with accretion 
rates > 10 -7 M /yr. In young dense cluster environments 
these episodal accretion processes will inevitably happen in 
addition to the steady viscous accretion of disc material. We 
showed that accretion-inducing encounters are more likely 
and lead to higher accretion rates for massive stars. 



Fig. 3. Average number of encounters in any time interval 
of 10 000 yrs that lead to accretion as a function of cluster 
age. 

believed to be powered by ep isodes of exceptionally rapi d 
accretion (up to 10 -4 M /yr) (|Hartmann fc Kenvon 1999) . 
Altern atively, such outbursts could be due to thermal insta- 
bility (|Bell fe Lin 1994t iBell 19991) . or to gr avitational in- 
stabilitv in the disc (jVorobvov fc Basu 2 006). While many 
of the properties (accretion rate, duration of burst) of the 
accretion bursts described here agree with those of FUors, 
the strongest argument against it is the lack of FUors 
observed in dense stellar regions. The arguments for and 
against a connection of encounter-triggered accretion bursts 
to FU Orionis objects will be discussed in more detail in 
Pfalzner (2008). 

Accretion bursts might be even more important in the 
earlier phases where star-disc systems are still in their for- 
mation process (jLarson 20021) . The simulation of the forma- 
tion of a small cluster of stars by Bate (2002) clearly shows 
the formation of many circumstellar discs. The complex 
dynamics of the system includes frequent protostellar in- 
teractions disrupting the discs by encounters. However new 
discs soon form from gas that continues to fall inwards, so 
that in this phase protostars may gain much of their mass 
in discrete episodes of encounter-triggered accretion. 

The observed accretion rates in Class I protostars 
are orders of magnitude smaller than those needed to 
form a star during th e lifetime of a Class I object 
(jMcKee fc Ostriker 20071 ). Kenyon et al. (1990) suggested 
two solutions to this problem - either significant accretion 
continues into the T Tauri stage or most of the accretion 
occurs in the embedded stage. The former appears to be 
ruled out because such stars accrete very slowly with no 
significant disc or envelope mass reservoir that they can 
draw on. In the latter case accretion could be episodic, so 
that the median accretion rate is much smaller than the 
mean. 

Another fact that can be explained if the accretion oc- 
curs in short bursts is that the observed luminosities of 
many protostars a re lower than is pred i cted for models wit h 
steady accretion (jKenvon et al. 19951 ICalvet et al. 2 004). 
The jet-like Herbig-Haro outflows probably powered by 
rapid accretion onto forming stars at early stages of evo- 
lution, are episodic or pul sed, suggesting that the accretion 
process is itself episodic ()Reipurth 200 1[ ). Remarkably the 
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Appendix A: Parameter study 

A.l. Encounter simulations 

In the encounter simulations themselves the mass of the 
discs was assumed to be m,j=0.001 M Q as is typical for T 
Tauri stars. However, the results are applicable for higher 
disc masses as well provided self-gravity effects do not play 
a dominant role. The main effect would be a rearranging 
of the disc material and a deviation from the Keplerian 
velocities of the disc particles. This probably happens for 
nid ~ M*. For such massive discs the cluster simulation 
particles would need to represent the mass of the star disc 
system and not just the stars as is currently the case. In 
this situation the encounter rate is likely to rise slightly in- 
creasing the number of encounter-induced accretion events 
even more. But this situation would need additional inves- 
tigations. 

In the encounter simulations the disc size was assumed 
to be 100 AU. But as the cluster consists of a wide spec- 
trum of stellar masses, the simulation results, valid for 
M-j* = 1 Mq , are generalized by scaling the disc radius ac- 
cording to rd = 15Qi/M*[M0], which is equivalent to the 
assumption of a fixed force at the disc boundary. The qual- 
ity of the dynamical models wer e determined by compar - 
ing them to observational data (|McCaughren et al. 2 002) 
at the approximate age of the ONC (l-2Myr). The quan- 
tities of interest were: number of stars, half-mass radius, 
number densities, velocity dispersion and projected den- 
sity profile. For mo re details of the selection process, see 
(|Olczak et al. 20061 ). 




Fig. A.l. Percentage of disc mass accreted onto star as 
function of the ratio of periastron r p to the disc radius rd 
for different M|/Mj* ratios in parabolic, prograde coplanar 
encounters. 
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4 


91 


4.03 


3.74 


3.14 


2.24 


1.83 


1.22 


0.54 


0.22 





1.0 


2.96 


5.01 


5.10 


4 


59 


4 


13 


3.17 


2.99 


2.23 


1.74 


1.61 


0.87 


0.30 


0.03 





1.5 


4.38 


4.51 


4.88 


4 


57 


2 


93 


2.35 


1.85 


1.46 


1.03 


0.87 


0.43 


0.01 








2.0 


4.97 


5.79 


4.69 


3 


36 


2 


09 


1.52 


1.38 


1.33 


0.74 


0.44 


0.06 











2.5 


5.76 


4.90 


4.03 


2 


40 


1 


69 


0.93 


0.79 


0.43 


0.17 


0.07 














3.0 


5.48 


4.11 


3.37 


1 


81 





98 


0.54 


0.47 


0.14 


0.03 


0.01 














3.5 


5.13 


3.60 


2.56 


1 


76 





63 


0.17 


0.04 


0.07 




















4.0 


5.37 


3.02 


1.91 





79 





14 


0.04 


0.01 




















0.0 


4.5 


4.60 


2.31 


1.53 





56 





11 





























5.0 


4.18 


1.76 


0.76 


8 


2 





























5.5 


4.09 


1.56 


0.73 


0.11 


































6.0 


3.75 


1.25 


0.38 


3 


































6.5 


3.69 


0.68 


0.27 







































7.0 


2.78 


0.45 


9 







































7.5 


2.64 


0.15 


3 







































8.0 


2.42 


9 


3 







































8.5 


2.11 


3 










































9.0 


1.60 


2 










































9.5 


1.24 


1 










































10. 


1.24 













































10.5 


0.55 













































11.0 


0.73 










0. 


































11.5 


0.36 













































12.0 


0.27 













































12.5 


0.17 













































13.0 


9 













































13.5 


8 


o 


o 












o 


o 


o 


o 


o 


(} 


o 


o 


o 


14.0 


7 













































14.5 


2 













































15.0 


2 













































15.5 
















































16.0 

















































Table A.l. Table of the percentage of particles accreted onto star 1 for all simulated configurations of prograde parabolic 
(e = 1) star-disc encounters. The first row contains the relative perturber masses M^/M^, the first column contains the 
relative periastra r p /r<i. 



